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Abstract. - We study a two-dimensional model for a long cylindrical stiff charged macroion 
immersed in a charge-asymmetric electrolyte with charge ratio +2/ — 1. The model is integrable 
and it allows an exact analytical determination of the effective charge of the macroion, which 
characterizes the electrostatic potential at large distances (compared to the screening length) 
from the macroion. At at high coulombic coupling, this model predicts charge inversion: for 
a highly negatively charged macroion, the effective charge could become positive, indicating 
an overscreening of the macroion by the divalent counterions. By studying the behavior of 
the coions and counterions density profiles close to the macroion, we show that the counterion 
condensation threshold is shifted to a lower value in absolute value. This plays an important 
role in the charge inversion phenomenon. 



The determination of the effective interactions between charged macroions immersed in 
electrolyte solutions is a central topic in colloidal science [1]. Based on the linear Debye- 
Hiickel (DH) theory, for a very long and thin (zero radius) cylindrical macroion, with linear 
charge density e/£, the reduced effective potential at a distance r is y(r) = 2XKo(Kr), where 
A = Ib/£, Kq is the modified Bessel function of order 0, y — etf)/(kBT) with tp the electric 
potential, e the elementary charge, T the temperature and ks is Boltzmann constant. Here, 
A is the reduced linear charge density, expressed in units of the inverse of the Bjerrum length 
Ib = e 2 / (ksTe), where e is the electric permittivity of the solvent. The Debye screening 
length is kT x = {AttIb Y^ a n a v a)^ 1 ^ 2 : where n h a is the bulk density of the microions of the 
electrolyte of type a and v a their valence. For highly charged macroions, the linear approach 
is inappropriate. A first improvement over the linear theory is to use the mean field nonlinear 
Poisson-Boltzmann (PB) equation. Under this approximation, the effective potential, at large 
distances, is again a screened potential: y(r) ~ 2X e gK (Kr) for nr ^> 1, but the prefactor A e fi 
is not anymore the bare charge A of the colloid, but it is known as the effective or renormalized 
charge [2,3]. 

PB approach is adequate [4,5] if the coulombic couplings between microions are small, for 
example for a two-component electrolyte: v\T <§C 1, vfr <§C 1, |wiU2|r <C 1 where T = 2ls/a, 
with a the average distance between microions in the electrolyte. For large coupling, new 
phenomena can occur that cannot be explained by the mean field PB approach. The most 
striking one is the phenomenon of charge inversion: a large fraction of counterions can condense 
into the macroion and the resulting dressed macroion can have a charge of opposite sign than 
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its bare charge. Closely related to this phenomenon is the possibility of attraction between two 
like-charged macroions. Several experimental observations and studies of these phenomena 
have been reported [6-8] and some theoretical explanations have been put forward [9-14], for 
a review see [15]. 

In this work, we report exact results for the renormalized charge of a cylindrical macroion, 
valid beyond the mean field approximation, in a whole range where the coulombic coupling can 
be large: < T < 2/3. The model under consideration is an infinitely long charged cylindrical 
colloid, with radius zero, immersed in a charge-asymmetric electrolyte solution composed of 
microions with valences V\ = 2 and v 2 = —1. Due to the translational symmetry along the 
direction of the colloid, we consider a two-dimensional (2D) model: we will be interested only 
in the microions density profiles and the electric potential along the radial direction from the 
colloid. Our model predicts the possibility of charge inversion: the effective charge and bare 
charge of the colloid have opposite signs when A < T — 1 < and T > 1/2. 

Being two-dimensional, our model describes correctly a cylindrical macroion immersed 
in an electrolyte formed by cylindrical parallel coions and counterions, with linear charge 
densities vi^e/a. As above, the coulombic coupling is T = 2/s/a. Many like-rod colloids and 
polyclcctrolytes can be described by the 2D model presented here. For example, synthetic 
polyelectrolytes such as poly(p-phenylene) and poly(styrenesulfonate) backbones used in water 
management have tunable linear charge density and a very large persistent length (20 nm) 
compared to other length scales involved in the problem. Of particular interest for biological 
and biophysics process, DNA (single and double stranded) and actin filaments are anionic 
cylindrical polymers that have also a large persistent length, and as such could be, in a first 
approach, described by the present model. With these cylindrical polyelectrolytes the situation 
described by this model could be reproduced experimentally. Experimental situations similar 
(although not exactly equal) to the one presented here are described in [7,16]. The experiments 
described in [7] and the possibility of charge inversion in those situations have important 
applications in gene therapy. Besides the possible applications to polyelectrolytes, it is well- 
known [17] that the 2D Coulomb gas is also a prototype model that can describe the physics 
of many interesting systems such as vortices in a superfluid, the XY model, and dislocations 
in a 2D crystal. 

For a cylindrical macroion with point-like coions and counterions, one important ingredient 
in the theoretical explanations [10] for charge inversion is the formation of a 2D liquid of coun- 
terions in the vicinity of the surface of the macroion and the strong longitudinal correlations 
between them. Our model, being two-dimensional, does not take into account these longitu- 
dinal correlations, but nevertheless shows that charge inversion is possible, thus showing that 
the radial correlations are enough to drive the charge inversion phenomenon. 

Our 2D model with zero radius charges is stable against the collapse of pairs of opposite 
charges provided that < T < 1. For technical reasons explained below our results are valid 
for < T < 2/3. To ensure the stability against collapse of counterions into the macroion, 
we require —1/2 < A < 1. The study of the effective charge beyond those limits requires 
to consider a macroion with nonzero radius and it is beyond the scope of this Letter, see 
however [18,19]. Provided —1/2 < A < 1, the introduction of a small hard core radius R (kR <C 
1) for the macroion is an irrelevant perturbation. These limits for A correspond to the Manning 
thresholds for counterion condensation [20] derived within a mean field approach. In all the 
theoretical approaches to charge inversion [9-11] the counterion condensation phenomenon 
is crucial. This is not an exception here. We will show that, due to the strong correlations 
between microions, the threshold for counterion condensation, in the case A < 0, is changed 
to (r — l)/2: it take place before than predicted by the traditional Manning theory. 

Exacts results for our model can be obtained due to the recent advances in the the- 
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ory of two-dimensional Coulomb systems: the exact bulk thermodynamics of the charge- 
symmetric [21] and the (+2/— 1) charge- asymmetric [22] two-dimensional two-component 
plasma are known. This has been possible due to their relationship with the integrable sine- 
Gordon (sG) and complex Bullough-Dodd (cBD) models. Using the field theoretical tools 
from the sG model, exact results for the effective charge and other quantities for a cylindrical 
macroion in a charge-symmetric electrolyte have been obtained [19]. Here, we use the known 
expressions for the form-factors of exponential fields of the cBD model to obtain results for 
the charge-asymmetric case. 

First, let us consider the electrolyte in the absence of any macroion. Carrying out a 
Hubbard-Stratonovich transformation, the grand canonical partition function of the elec- 
trolyte can be cast as the partition function of the cBD model, with action [22] 
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where Zi t 2 are the fugacities of the microions and b= VT. To give a precise meaning to the 
fugacities in the cBD model the conformal normalization should be used: (</>(0)<^(r)) Zl=Z2=0 = 
— lnr. In the following we shall denote (• • •) an average with the cBD action (1). We now 
consider that a single macroion with charge A is immersed in the electrolyte at the origin. Let 
Q = A/T. The density of positive and negative microions at a position r from the macroion 
arc [19,23,24] 

/ p ibQ<j>{0) p ibq±<t>(r)\ 
n ± ( r ) = n ± , TnZT! h ^ ( 2 ) 

v ' ± lp}bQ(j>\ l e ibq±<p\ V ' 

with q + = 1 and g_ = —1/2 respectively and n b ± are the bulk densities, far from the colloid. 
The above correlation function can be expressed as a sum over all intermediate iV-particle 
states of the cBD model as [22,25] 
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where the labels the particle spectrum of the cBD model. The rapidity of the i-th particle, 
which is of the type e,, is 6*, and has mass m ei . Fq and F q± are the form- factors of exponential 
fields in an A^-particle state with particle spectrum ei, . . . , ejv. The particle spectrum of the 
cBD is studied in [26]. This expression is appropriate to find the large distance expansion of 
the correlation functions. The dominant term is obtained by considering the lightest particle 
in the spectrum of the cBD model. For T small enough (see below) the lightest particle is the 
1-breather. Its mass is known [27] and it can be expressed in terms of the Debye length as [22] 
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where £ = T/(2 — V). The corresponding form-factor is [28,29] 
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and 

j = _ 4 r cosh(f)sinh(f)sinh(^) rft 

fc 7 sinhtcosh(i/2) i [ ' 

Notice that the form- factor diverges for T = 2/3. This indicates a change of behavior in 
the expansion (3): some subdominant terms (for T < 2/3) become of the same order as the 
contribution of the 1-breather at T — 2/3 (see [24] for a similar situation in the short-distance 
expansion of correlation functions of the sG model) . In the following, we restrict our analysis 
to T < 2/3. 

Let us define the effective interaction energy (in units of ksT) E\^ q± (r) of a microion 
with charge 2q± (+2, —1 respectively) of the electrolyte with the macroion by the relation 
n±(r) = n b ± exp[— E\x q± (?*)]. At large distances, r — > oo, n±(r) ~ n±(l — E\.r q± (r)). Using 
the expansion (3) with only the dominant contribution from the 1-breather form-factor, we 
obtain, for mr 3> 1, 

8v/3cos(f (1 - 20) h ^ ( 2ir\ \ _ f 2n\ tt tt£ 



~ — ^f(TT# ^ sin ^3(2-T) J cos IsT^T) - 6 " f ) Ko{mr) (8) 

8V3 cos (f( 1 + 4 ^)) sin (¥) r . ( 2tt\ \ ( 2tt\ tt ttA Tjr , 
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(9) 

Replacing these expressions into the density profiles and integrating Poisson equation, we find 
the effective electrostatic potential y(r) created at a distance r S> vnT 1 from the macroion. 
At large distances, the electrostatic potential takes a similar form to the one predicted by 
the linear DH theory y(r) ~ 2X c ffK n (mr), supporting the hypothesis of the existence of an 
effective or rcnormalized charge A c ff. Notice however that there is also a "renormalization" 
of the screening length: it is mT 1 given in (4), instead of the usual Debye length For 
< T < 2/3, the screening length m" 1 is smaller than the Debye length The effective 

charge reads 

CV3 S in(7rQe J > sin (jg^) cos (jg^ - g - f ) ^ 



3 sin (f (1 + 0) sin sin (f (1 + 0) sin 



In the low coupling limit r « 1 we recover the results from PB theory [18,30]: A c ff = 
\/3[2sin(^ - f ) + 1]/(2tt) as expected. 

Figure 1 shows the effective charge as a function of the bare charge. An important result 
is that the phenomenon of charge inversion is possible: when T > 1/2 the effective charge 
and the bare charge have opposite signs when — l/2<A<r — 1<0. The charge inversion 
occurs for a highly charged negative macroion which is overscreened by positive microions of 
valence +2 resulting in a positive effective charge. On the other hand, for a positive macroion, 
screened by microions of valence —1 the charge inversion does not occurs. 

The mean field PB formalism does not predict the charge inversion phenomenon, thus 
the microions correlations are fundamental for this phenomenon to take place. The charge 
inversion has previously been predicted for a planar charged interface [12] and for spherical 
colloids [13] with divalent counterions and monovalent coions within the hyperncttcd chain 
approximation (HNC). Our model, which is exact (we have solved exactly, without any ap- 
proximations, the statistical mechanics of the model), put on more firm ground the predictions 
of [12, 13] and thus suggests that the HNC approximation captures adequately the microions 
correlations responsible of the charge inversion. 
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Fig. 1 - The effective charge A c ff as a function of the bare charge A of the colloid, for different values 
of the coupling constant T and for the PB approximation. The inset details the cases V — 0.5 and 
r = 0.6 for A < 0. This later case shows the charge inversion phenomenon. 

The charge inversion can be actually very large: for example, for T = 0.65 and A = —0.48 
we have A e ff = 1.22, thus a charge inversion ratio of more than 200% . This is similar to the 
"giant" charge inversion studied in [31], although here the effective charge characterizes the 
long distance signature of the potential, contrary to its short-distance behavior as considered 
in [31]. The charge asymmetry of the electrolyte seems to be fundamental for the charge 
inversion to take place. Exact results for a symmetric electrolyte show no charge inversion [19, 



The charge inversion is accompanied by a change of behavior in the large distance behavior 
of the effective interaction between the macroion and the microions. For A > T — 1, at 
large distances, the effective interaction between the macroion and a counterion is attractive, 
whereas its interaction with a coion is repulsive. However, as it can be seen from (8-9), when 
r > 1/2 and — l/2<A<r — 1<0 the interaction of the macroion with a counterion is now 
repulsive and its interaction with a coion is attractive, provided mr 3> 1. 

A similar change of behavior also occurs at short distances, but only for the interaction 
with a coion. We can obtain the short-distance behavior of the density profiles by using the 
operator product expansion (OPE) in the correlation function appearing in (2). The OPE 
for the cBD model has been developed in [32]. Using the OPE (see also [24,33]) we find, 
for — 1/2 < A < that n+(r) oc r 4A when r — > 0, yielding an attractive effective interaction 
Ex,+r(r) ~ — 4Alnr of the macroion with a counterion (charge +2) at short distances, which 
is the expected behavior. On the other hand the density profile of the coions (charge —1) at 
short distances exhibits a change of behavior: 



23]. 



n_(r) oc 




for S^i < A < 

for - I<A<^-<0 



(11) 



giving an effective potential at short distances, mr <C 1 
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The change of behavior at A = (r — l)/2 can be understood as a first step in the counterion 
condensation [23,24] for large coupling T. A fraction of the counterions are condensed into the 
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macroion so that the effective interaction of the macroion with the coions at short distances is 
not anymore the bare Coulomb potential 2Alnr, but it changes to 2 (T — 1 — A) lnr. The strong 
correlations between the microions shift the threshold for counterion condensation from the 
Manning value — 1/2 to (r — l)/2. Since the trademark of the condensation is first noticed in 
the change of behavior (11) of the coions density profile, this shift happens only for a situation 
with added salt (contrary to the no salt case [34]). When r — > 0, the above limit reduces to 
the standard value A = —1/2 [20]. In PB theory, the effective interaction between a coion 
and the macroion behave et short-distances as the bare Coulomb potential 2Alnr down to 
A = -1/2 [18]. 

Interestingly, as A decreases beyond A < (r — l)/2 the effective interaction, at short 
distances, of the macroion with the coions becomes less and less repulsive. Paradoxically, it 
even becomes attractive when A < T — 1, provided T > 1/2. The divalent counterions are 
strongly attracted to the macroion and since they attract the monovalent coions, there can be 
a net attraction of the coions to the macroion. This is accompanied, as we have seen before, 
with the charge inversion at large distances. The very different behavior, at high coupling, of 
the coion and counterion density profiles at short distances shows that it is difficult to define a 
"dressed" charge of the macroion which characterizes the potential at short distances. Indeed, 
when A < (T—l)/2 < 0, the coions see a dressed charge T — 1 — A which can be even positive, 
whereas the counterions still see a charge A < 0. 

Returning to the charge inversion phenomenon, we should point out that it does not 
necessarily implies attraction between like-charge macroions. In particular for two identical 
macroions with charge A < T — 1<0, as both macroions have an inverted effective charge, 
the net effect would still be a repulsion between the macroions. To study the possibility of 
like-charge attraction, let us consider two macroions with linear charge densities Ai and A 2 
immersed in the electrolyte at the origin and at r respectively. In the language of the cBD 
model, the effective interaction between the two macroions is given by [19,23,24] 

/ e ibQ 1 4>(0) e ibQ 2 4>(r)\ 

exp [-E XlXa (r)] = {elbQlfP){eibQ2<P) (13) 
with <3i,2 = Ai, 2 /r. Using the form-factor theory explained above we find, for mr 3> 1, [22] 

with 

r . / 27rAi, 2 \ / 27tA 1>2 tt ttA 

Al - 2 = sl H372^J cos l372^T)-6-yJ • (15) 

The "charges" Ai i2 exhibit the same changes of sign as the effective charge (10) shown in 
Fig. 1. Thus, two negatively charged macroions can have an attractive effective interaction at 
large distances provided that Ai < T — 1 < and r — 1 < A 2 < for T > 1/2. Also a positive 
macroion (say A 2 > 0) can have a repulsive effective interaction with a negative macroion 
(charge Ai) provided Ai < T — 1 < 0. 

Summarizing, we have considered a model that describes the effective charge of cylindrical 
macroions immersed in a charge- asymmetric +2/— 1 electrolyte solution. This model is ex- 
actly solvable and predicts charge inversion at high coupling T > 1/2 for negatively charged 
macroions provided their linear charge density A < T — 1 < 0. Also, we have shown that a shift 
in the counterion condensation threshold from A = —1/2 to A = (r — l)/2 occurs due to the 
strong coupling between the microions, which is part of the mechanism of charge inversion. 
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